Optimized sympathetic cooling of atomic mixtures via fast adiabatic strategies 
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We discuss fast frictionless cooling techniques in the framework of sympathetic cooling of cold 
atomic mixtures. It is argued that optimal cooling of an atomic species - in which the deep- 
est quantum degeneracy regime is achieved - may be obtained by means of sympathetic cooling 
with another species whose trapping frequency is dynamically changed to maintain constancy of 
the Lewis- Riesenfeld adiabatic invariant. Advantages and limitations of this cooling strategy are 
discussed, with particular regard to the possibility of cooling Fermi gases to a deeper degenerate 
regime. 
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The continuous progress in achieving low temperatures 
has opened up new research directions in atomic physics, 
most notably the formation and study of degenerate Bose 
and Fermi gases. For the latter, most of the interesting 
physics, such as novel supcrfluid states or ordered phases, 
is conditional on the achievement of very low degeneracy 
factors T/Tp, now limited to above ~ 1CT 2 [10. In 
most of the experiments the Fermi gas is in thermal con- 
tact with a Bose gas cooled with the usual evaporative 
cooling strategy and having a much larger heat capacity 
than the Fermi gas, resulting in the so-called sympathetic 
cooling (see [3|-|5[ for earlier demonstrations). However, 
both the loss of atoms intrinsic to evaporative cooling 
and the strong dependence on temperature of the heat 
capacity limit the efficiency of sympathetic cooling 0, Q . 

If the main goal of the experiments is to achieve the 
deepest Fermi degeneracy, an alternative procedure can 
be envisaged with the use of species-selective traps in 
which the Bose gas alone is kept as close as possible to 
the nondegenerate regime, as proposed for bichromatic 
optical traps in 0], and for light-assisted magnetic traps 
in Q. Recently, species-selective light-assisted magnetic 
traps have been successfully implemented to sympathet- 
ically cool Bose-Bose mixtures such as 41 K- 87 Rb and 
87 Rb- 174 Yb [llj, and these experimental techniques could 
be naturally extended to Fermi-Bose mixtures. However, 
the ultimate, optimal ideal coolant would be a gas that 
retains the maximum value of its heat capacity, both be- 
cause it does not enter the degenerate regime, and also 
because the number of atoms is not diminished during 
the cooling stage. These two features are common to a 
new class of cooling techniques known under the name 
of frictionless (or, less rigorously - as seen in later con- 
siderations - fast adiabatic) cooling. This method was 
pioneered in atomic physics in fllj ]. and experimentally 
implemented for fast decompression of 87 Rb atoms both 
in the nondegenerate [Hj] and degenerate flp regimes, 
and for demonstrating fast atomic transport [14[. 

It is the purpose of this Rapid Communication to ex- 



plore the extent to which frictionless cooling techniques 
may be useful in sympathetic cooling of Fermi gases, 
while also discussing their practical limitations. An im- 
portant motivating factor is that an usually undesired 
feature of frictionless cooling, i.e. the fact that the 
atomic cloud does not increase its phase space density 
and therefore its degeneracy, is instead a crucial asset 
from the perspective of maintaining the gas in the non- 
degenerate regime and making it an ideal coolant. 

Fast frictionless cooling ll| is based on the use of 



an adiabatic invariant for harmonically trapped atoms 
of mass m - hereafter assumed in the nondegenerate 
regime and therefore with negligible interactions. Fol- 
lowing Lewis and Riesenfeld [la, [la] , we consider the op- 
erator I(t) = t: 2 /2m + mui 2 q 2 / (2b 2 ) where 7r = bp — mbq 
is the momentum operator conjugate to the operator 
q/b. The parameter loq can be chosen as the initial fre- 
quency, and b(t) is a time-dependent frequency scaling 
factor which must satisfy, in order for I{t) to be an in- 
variant operator, the Ermakov equation b(t) + ui(t) 2 b = 
Lu 2 /b 3 . This can be solved by imposing boundary con- 
ditions on b(t) and its first and second time derivatives, 
if we assume both a targeted final trapping frequency 
Wf and a total time duration for the adiabatic protocol 
U U( , thus obtaining a well-defined Ermakov trajectory 
for the requested time-dependent trapping frequency as 
uj\t)=^ Q /b\t)-b{t)/b{t). 

In our analysis, we compare three decompression 
strategies for an harmonically trapped atomic cloud as 
shown in Fig. 1 where, for concretcness, we assume 
the initial trap frequency loq — 2it x 250 Hz and the 
final angular frequency u>{ = 2ir x 2.5 Hz, ujq/lu{ = 10 2 . 
Two of the trajectories for oj 2 (t) (continuous and dashed 
lines) are solutions of the Ermakov equation and there- 
fore leave the Lewis- Riesenfeld operator invariant in time. 
The main difference in the two trajectories is the absence 
(continuous) or the presence (dashed) of a stage in which 



the square frequency becomes negative. As shown in [11 



the presence of an anti-trapping stage appears for small 
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FIG. 1: (Color online) Comparison of various decompression 
strategies as seen from the square of the trapping frequency, 
normalized to its initial value, versus time normalized to the 
overall duration of the decompression. The curves here and 
in the following figures indicate an Ermakov trajectory with 
positive square frequency lasting 25 ms (black, continuous 
line), a fast Ermakov trajectory of total duration of 6 ms 
(red, dashed line), and a linear ramp down of the trapping 
frequency approximating a conventional slow adiabatic strat- 
egy (blue, dot-dashed line) with duration of 400 ms. Notice 
in particular, that ui(t) 2 /ujq < in a finite interval in the 6 
ms case, leading to an anti-trapping stage, which occurs for 
all choices of tf < 25 ms. 

targeted times tf. Physically, the anti-trapping stage en- 
ables the wave function to spread out faster so as to reach 
the final target width within given tf. It is interesting 
to note that the Ermakov equation typically produces a 
trajectory where, after anti-trapping, the trap frequency 
first increases to a value above Uf before finally relaxing 
down to ojf. This may be interpreted as a way to con- 
trol the excess momentum directed away from the origin 
attained during the anti-trapping stage. Experimentally, 
an anti-trapping stage can be obtained by properly mod- 
ulating the power of blue-detuned beams, as discussed 
in Q. For comparison, the third trajectory shown (dot- 
dashed line) is a naive decompression strategy based on 
a linear ramp down of the frequency to its targeted fi- 
nal value. Its overall duration is chosen as tf = T max , 
where T max = 2-n/uif is the longest harmonic oscillator 
period during the evolution. In this way we compare two 
Ermakov-rclatcd cooling strategies to a strategy which 
is adiabatic in the usual sense of the word, i.e. with 
parameters changing slowly with respect to the intrinsic 
timcscale set by the harmonic oscillator frequencies. 

We have numerically solved the Schrodinger equation 
for the time-dependent harmonic oscillator evolving un- 
der the u(t) solution of the Ermakov equation. By com- 
paring the solution with the solution obtained from the 
effective Gaussian dynamics formalism [ItJ it was found 
that the wave function, if prepared initially as a Gaus- 
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FIG. 2: (Color online) Thermodynamics of the decompres- 
sion. Temperature scaled to initial temperature (top) and 
Shannon entropy (bottom) versus time for the three cooling 
strategies described in Fig. 1. In the case of the decompres- 
sion strategy corresponding to tf — 6 ms the temperature 
exceeds the initial temperature in a time interval related with 
the time interval in which the square frequency becomes neg- 
ative, while in the case of tf = 25 ms the temperature be- 
havior is decreasing monotonically. The threshold situation 
in which the maximum temperature during the anti-trapping 
stage equals the initial temperature is achieved for tf « 11 
ms. 



sian, retains its Gaussian form throughout the evolution, 
which proves useful in characterizing the temperature of 
the cloud. The velocity distribution is routinely used 
to estimate experimentally the temperature of ultracold 
atomic gases. In our case, the Gaussian state can be 
used to represent an ensemble of atoms with a normally 
distributed position and momentum such that, by as- 
suming Maxwell-Boltzmann distribution, the variance of 
the velocity distribution is er^ = fc^T/m where fc^ is 
the Boltzmann constant and T is the temperature of the 
atomic cloud. At the same time, the use of the variance 
of the velocity distribution, or equivalently the momen- 
tum variance, allows us to monitor the energy population 
probabilities during the cooling procedure. This is simply 
a reflection of the fact that a larger momentum variance 
implies the occupation of a large number of excited mo- 
mentum states. This complements Ref. [ll[ where the 
total energy is used as an indicator of temperature. The 
time evolution of the temperature is shown for the three 
cases in the top plot of Fig. 2. In addition, the en- 
tropy S(t) is evaluated based on its Shannon definition, 
S(t) = — J2 n \ c n(t)\ 2 l°g l c «W| 2 where the instantaneous 
cigenstate basis \ijjn{t)) corresponding to the harmonic 
trap frequency uj(t) is used in the projection of the time- 
evolved wave function. Although S(t) is the Shannon 



entropy instead of thermodynamic entropy, an adiabatic 
process would keep c„(i) constant in time, meaning the 
entropy would remain constant. 

In the case of the frictionless strategy with an anti- 
trapping stage the basis for the inverted harmonic oscil- 
lator [18| was used (h = m = 1): 
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where tp± (x, t) corresponds to the solution for which the 
particle is confined to an expanding box with a moving 
boundary condition. TV is the normalization constant and 



the discrete energy eigenvalues are E n 



e where 



e = (?i + i) 2 7r 2 /Lg with n an integer and L is the width 
of the wave function as it first enters the inverted po- 
tential. ip2(x,t) corresponds to the plane wave solution 
with k being the wave number, and since the plane wave 
solution cannot be normalized, there is no normalization 
constant. It was found that, as the anti-trapping stage 
lasts for extremely short duration typically around 0.25if , 
there is not enough time for our wave functions to cou- 
ple to the free- particle plane wave solution of tjj2(x,t). 
Therefore our wave functions project entirely onto eigen- 
states ipi (x, i). The results of this entropy calculation are 
shown in the bottom panel of Fig. 2. 

Two general features are worth commenting on. First, 
for the two cases corresponding to the conservation of 
the Lewis- Riesenfeld invariant, entropy is clearly not con- 
stant, which implies that adiabaticity does not hold dur- 
ing the cooling process. The maximum entropy appears 
at the point where the trap is shallowest while in cases 
where an anti-trapping region exists, the spikes in en- 
tropy are where the curvature of the trap changes its sign. 
However, the final entropy exactly equals the initial one 
for these cases, as if a real adiabatic process has taken 
place. This is in stark contrast to the linear ramp down 
decompression where the entropy is constant throughout 
the evolution except as we approach the final time. This 
increase is expected as the decreasing frequency leads to 
stricter conditions for true adiabaticity, indicating a need 
for much longer timescales for the ramp-down. 

Second, the difference between the linear decrease and 
the Ermakov constructs is also clear in the temperature 
dynamics. Higher intermediate temperatures correspond 
to the presence of an anti-trapping stage, which in turn 
appears in the Ermakov solution for shorter if. We note 
that in general the presence of a stage in which the tem- 
perature of the cloud may become much higher than the 
initial cloud temperature is problematic. Since atomic 
traps have a finite potential depth, a rise of the cloud tem- 
perature will induce losses of the most energetic atoms 
UJ. From our standpoint, this leads to a practical limita- 
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FIG. 3: (Color online) Dependence of the maximum tem- 
perature attained versus the duration of the cooling strategy 
tf for the targeted temperature decrease T{tt) /T(0) = 10~ 2 . 
Solid line: maximum achieved temperature scaled to the ini- 
tial temperature. Dashed line: maximum temperature ex- 
pressed in units of quanta of the harmonic oscillator with the 
instantaneous frequency corresponding to the time at which 
the maximum temperature occurs. For tf larger than ~ 11 
ms, the initial temperature is also the maximum temperature. 
For shorter tf, substantially higher temperatures occur after 
the anti-trapping region as seen in the top plot in Fig. 2. 



tion on the minimum possible time for the cooling strat- 
egy, with a natural tradeoff also determined by the life- 
time of the atoms in the trap. The final temperature for 
the linear ramp case is slightly higher than that for the 
other two cases. However, if the linear ramp occurred on 
the same time scales as the other two fast adiabatic meth- 
ods, the final temperature would be much higher [20| . In 
order to quantify this potential bottleneck of the cool- 
ing protocol, we show in Fig. 3 the maximum achieved 
temperature, i.e. the peak value of the temperature in 
Fig. 2, scaled to the initial temperature, versus the time 
duration of the cooling procedure. In the same figure we 
also show the maximum achieved temperature in units 
of quanta of the corresponding instantaneous harmonic 
oscillator ksT max / huj(t meix ) . This figure of merit is use- 
ful to assess trap losses in the case of a finite trap depth. 
The results show, that for tf shorter than approximately 
11 ms, the maximum temperature of the atoms becomes 
comparable to the energy of the 15 th excited state of the 
instantaneous harmonic trap, and quickly increases for 
shorter tf. Nevertheless, it appears that optimal cooling 
with minimal trap loss on timescales comparable to the 
oscillation period of the trap is possible as long as exces- 
sively short cooling times are avoided. This is clear from 
the temperature curve where for Tf > 11 ms, the initial 
temperature is also the highest temperature indicating 
stability in the dynamics as the trap is made more shal- 
low. By contrast, at shorter times, substantial heating re- 
sulting from the anti-trapping region would make the fric- 
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tionless strategy difficult to implement. The maximum 
heating due to anti-trapping was found to occur typically 
at t ~ if/2, the time when the trap changes its curvature 
at the end of anti-trapping. The sudden "jump" occurs at 
around if sa 11 ms, since for this case the maximum heat- 
ing due to anti-trapping has the same magnitude as the 
initial temperature. We emphasize that, as anti-trapping 
occurs for all tf < 25 ms in our example, the presence of 
an anti-trapping stage alone is not a sufficient condition 
for heating of the atoms to temperatures greater than the 
initial temperature. 

A different issue arises when the spatial overlap of 
the two clouds is taken into account, as the coolant 
species considerably increases its size during its decom- 
As discussed in [11] . the position variance 



prcssion. 

is directly related to the frequency scaling factor b as 
a\ = h(n + l/2)t 2 /(mw ). For the two adiabatic in- 
variant strategies, the temporal variation of the position 
variance is independent of t[ as expected from the fact 
that the position variance is proportional to b(t). Since 

^(*f)/°"x(°) = fo2 (if)/ &2 (0) = WoM = 10 2 as a result of 
the boundary condition for b(t) , it is conceivable that if a 
large b(t{) is targeted, the spread of the atomic cloud of 
the coolant will result in a small overlap with the cloud to 
be sympathetically cooled. The linear ramp down result 
also exhibits a large broadening expected of an adiabatic 
process that relaxes the trap frequency by a factor of 
10 2 . The overlap issue seems to be crucial in various 
experiments (see for instance [2l[ for the effect of grav- 
itational sagging on a large mass ratio Fermi-Bose mix- 
ture), and its effect on species-selective traps has been 



discussed in detail in |22j. We note that, in practice, 



the issue of spatial overlap is less severe than expected, 
since for accidental reasons the fermionic isotope of the 
alkali primarily used in the experiments, 6 Li, is lighter 
than the bosonic species used as coolants, 23 Na, 87 Rb, 
and bosonic isotopes of Yb. Therefore the initial over- 
lap between the two species sees the fermionic species 
more spread out than the bosonic species. The decom- 
pression of the bosonic species alone in general allows for 
a better spatial overlap at intermediate times, but this 
also provides a limit on the maximum allowed decompres- 
sion before the spatial overlap decreases again. It is also 
worth to point out that, especially for very short cool- 
ing times, a possible issue arises with the sympathetic 
equilibration rate between the two species, depending on 
the interspecies elastic scattering rate. Although this is- 
sue should be studied in detail for each specific atomic 
mixture, in general this is less relevant than the spa- 
tial overlap, and can be circumvented by using magnetic 
or optical Feshbach resonances to boost the interspecies 
elastic scattering length as in the case of the 6 Li- 87 Rb 
Fermi-Bose mixture [23l - l25l j. 

In conclusion, we have discussed the possibility of 
achieving deep Fermi degeneracy via sympathetic cool- 
ing using a frictionless strategy for the coolant species. 



The identified advantages in this setting are the maximal 
heat capacity retained by the coolant due to the conser- 
vation of the number of atoms and the preservation of its 
phase space density in the nondegencrate regime where 
the specific heat retains its Dulong-Petit value. At the 
same time we have also identified two important limi- 
tations, a cooling stage in which the temperature rises 
significantly allowing for trap losses if the potential en- 
ergy depth is not large enough, and the spreading of the 
cooling cloud which will reduce the spatial overlap with 
the fermionic cloud. The coolant does not need necessar- 
ily to be a Bose gas, as it is kept in the non-degenerate 
regime during the entire process, thus providing further 
flexibility in any possible experimental implementation 
of this scheme. 
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